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P and CP violation in cosmology can be manifested as large-scale helical velocity flows in the ambi¬ 
ent plasma and as primordial helical magnetic helds. We show that kinetic helicity at last scattering 
leads to temperature-polarization correlations {Cf^ and Cf^) in the cosmic microwave background 
radiation (CMBR) and calculate the magnitude of the effect. Helical primordial magnetic helds, 
expected from cosmic events such as electroweak baryogenesis, can lead to helical velocity flows and 
hence to non-vanishing correlations of the temperature and B-type polarization. However we show 
that the magnitude of the induced helical how is unobservably small because the helical component 
of a magnetic held is almost force-free. We discuss an alternate scheme for extracting the helicity 
of a stochastically homogeneous and isotropic primordial magnetic held using observations of the 
CMBR. The scheme involves constructing Faraday rotation measure maps of the CMBR and thus 
determining the sum of the helical and non-helical components of the primordial magnetic held. 
The power spectrum of B-type polarization huctuations, on the other hand, are sensitive only to 
the non-helical component of the primordial magnetic held. The primordial magnetic helicity can 
then be derived by combining these two sets of observations. 

PACS numbers: 98.80.Cq 


The cosmic microwave background radiation (CMBR) 
has become a remarkable tool for mapping the past and 
present states of the universe. As more refined observa¬ 
tions are made, a larger array of theoretical ideas will 
be put to the test and more details of the history of the 
universe will emerge. 

Until now, only the temperature anisotropy of the 
CMBR has been measured. It is widely recognized that 
even the most precise measurement of the temperature 
spectrum would still be compatible with a large class 
of initial conditions for the fluctuations in the primor¬ 
dial universe. Measurements of the CMBR polariza¬ 
tion would signihcantly reduce this degeneracy. The 
much anticipated hrst detection of the CMBR polariza¬ 
tion could be provided by the MAP satellite, launched by 
NASA in June of 2001, or balloon-borne B2K and MAX- 
IPOL experiments. There are several other ongoing or 
planned ground-based experiments that have the detec¬ 
tion of CMB polarization as their primary goal. How¬ 
ever, most expectations of obtaining polarization spectra 
are connected with the Planck satellite mission, currently 
scheduled to be launched in 2007. 

A possibility that has already received some attention 
is that large-scale Parity (P) violation may be observed 
via the CMBR ^ . In Ref. Q it was shown that a co¬ 
herent magnetic field would induce non-zero P-violating 
correlations in the CMBR through Faraday rotation. In 
Ref. (H the P violation was due to the dynamics of a pseu¬ 


doscalar field which could induce a fixed-handed rotation 
of the polarization plane of the CMBR. In another case 
studied in Ref. the pseudoscalar held could produce 
an excess of left- over right-handed polarization states of 
gravitational waves, which would be manifested in the 
CMBR as a P-violating signature. In Sec. I we exam¬ 
ine the consequence of yet another possible source of P- 
violation (and also CP-violation), namely large-scale heli¬ 
cal velocity hows {i.e. kinetic helicity) at recombination. 
We hnd that the CMBR temperature and polarization 
cross-correlators can provide a measurement of kinetic 
helicity (Sec. ||). A motivation for considering helical 
hows is that helical magnetic helds mw have been gen¬ 
erated in the early universe ill I and our initial 
guess was that such helds could induce helical velocities. 
However, in Sec. HI we show that the magnitude of the 
kinetic helicity induced by helical magnetic helds is in- 
signihcant and does not lead to any observable signature 
in the CMBR. 

This leaves us with the challenge of devising a method¬ 
ology to observe the helicity of primordial magnetic helds. 
In Sec. 0 we describe a method which involves observ¬ 
ing the temperature and polarization huctuations of the 
CMBR and also the Faraday rotation measure of the 
CMBR polarization along different directions. From a 
suitable combination of these observations, it is possible 
to extract the primordial magnetic held helicity. Hence, 
at least in principle, the helicity of cosmic magnetic helds 
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may be observable. 


I. KINETIC HELICITY 


Py (k) to be a tensorial function of fcb The latter con¬ 
dition is needed to ensure that the correlator Pij(k) de¬ 
scribes a homogeneous and isotropic random vector field. 
It follows from these restrictions that 


Kinetic helicity of a velocity field v is characterized 
by a non-vanishing value of (v • a;), where u; = V x 
V. For the purpose of calculating the effect of helical 
flows on the cosmic microwave background radiation, it 
is useful to decompose the velocity field at last scattering 
into gradient components and rotational components 
as follows: 


Vj (x) 


f (fk 

~ J 


where 


(k) 


ei ± ie2 
-F— 

V2 


( 1 ) 

( 2 ) 


and the vectors ei 2 are such functions of k that they 
form an orthonormal basis together with 63 f k. The 
functions can be thought of as eigenvectors of he¬ 
licity because ik x where s = ±1. Also, 

k • = 0 and (k) = Q^"") (-k). Under a par¬ 

ity transformation, v(x) ^ —v(—x), k ^ —k while 
and therefore (k) ^ (—k), 

.j;(±l) —j. ( —k). 

In Appendix 1^ we will show that the parity-odd CMB 
correlators C^^^and are linearly dependent on the 
expectation value of the parity-odd combination of : 

(k) (-k) - (k) (-k)) . (3) 


The average helicity of the velocity field is proportional 
to the same quadratic combination of Namely, 

I d^x v(x + y) • [V X v(x)] = I 
X (u(i)(k)u(i)(-k) - z;(-i)(k)u(-i)(-k)) , ( 4 ) 

where k = |k|. 

Let us assume that the initial velocity field is random 
with a given power spectrum, such as 

k^ 

{v.{k)v,{k'))=v^,^{2nf6{k + k')P,,{k), (5) 

where vq is the characteristic velocity of the flow and 
fc* is the wave vector corresponding to a cutoff scale. 
In general, Py does not have to be a symmetric ten¬ 
sor. The restrictions on it are: the reality condition 
Pj* (k) = Py(—k) = Pji{k), the divergence-less condi¬ 
tion (incompressibility) k^Pij = 0 , and the requirement 
of correct transformation under rotations, which forces 


Pvi^) = fik) 


Sij ki kj 


■ ig{k)sijik‘, 


( 6 ) 


where f{k) and g{k) are real functions of fc = |k|. By 
using 


(|a;(k)±A:v(k)|2)>0, 

where w(k) is the Fourier transform of <^(x), we can de¬ 
duce that the functions / and g are not completely inde¬ 
pendent but must satisfy the inequality 

/(k)>|5 (k)|. (7) 

[/(k) > 0 follows by taking the trace of Eq. (||).] Eqs. (|^) 
and i) now lead to 

(k) (k') — (k) (k^^ 

= u2^(27r)3<5(k + k')P*^(k)Q[,'^(k)0;f'^(k). (8) 

/vjft 

The /-term in Eq. (jbj) will not produce helicity because 
it is symmetric in i, j. Since 

ie,,ik^Qf{k)Q^-^\k) = l, (9) 

we get 

un 

= vl^{2^)H{k + \^)g{k). ( 10 ) 

/C;)5 


II. CMBR SIGNATURE OF KINETIC 
HELICITY 

The CMBR temperature and polarization at a given 
direction on the sky are described by time-averaged com¬ 
ponents of the electric field intensity tensor lij = EiEj. 
Equivalently, the CMBR is represented by the Stokes pa¬ 
rameters I, Q, U and U |^, ||, where I is the total inten¬ 
sity, Q and U are two components of the linear polariza¬ 
tion and V quantifies the circular polarization. Assuming 
that the CMB photons prior to the last scattering were 
unpolarized, we can drop V from consideration, since 
Thomson scattering can only generate linear polariza¬ 
tion. 

Decomposing the temperature anisotropy into spheri¬ 
cal harmonics gives 

A (n) = '^ ^ ^ ^ Yi^ (n) af^ . (11) 
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FIG. 1: {1/2tx)1{1 + l)C'P’-®®(fc.»7o)"+V^^o vs Z for n = 2 

and 271 jk, = 288 Mpc, where n, fc* and vo are defined in 
Eq. §) and 770 is the conformal time today. We have also 
taken g(k) = 1 [defined in Eq. (^]. The solid line shows the 
TB correlation and the dashed line is EB. Dotted lines show 
the same spectra in the absence of a cutoff. 


Analogously, Q (fi) and U (n) can be decomposed using 
spin-2 spherical harmonics (see Refs. IJ): 

= ( 12 ) 
l,m 

The “electric” and “magnetic” components of polariza¬ 
tion are eigenstates of parity and may be defined by 

= -i (-2)\ B ^ _i_ /„(2) 

'^Im 2 \^lm ' ^Im J ’ ^Im 27 \^lm ^Im J ' 

(13) 

Under parity inversion the components transform as 
“im ^ (-1) Oim . while ai^ - (-1) 

Observations of CMBR are usually presented in the 
form of spectral functions Ci ’ defined by 

E (14) 

m— — l 

where X and Y stand for T, E or B. Correlators in 
Eq. ( 0 ) are real because they involve summation over 
both positive and negative values of m. This can be 
seen by noting that the tensor spherical harmonics sYim 
satisfy 

[.Yirr,]* = (-1)*+™ _sYi^_^ , (15) 

and therefore the components a^, 0,'h^’^ transform un¬ 


der complex conjugation as 




= (-ir 


a 


=F2 

l,—m ’ 


T,E,B 
‘I 


(-ir 


a 


T,E,B 
l, — m ■ 


(16) 

The correlators Cf^ and are parity-odd, while 

all other correlators are parity-even. If no parity- 

odd sources are present, the ensemble of fluctuations is 
statistically parity-symmetric and therefore the correla¬ 
tors and must vanish. This is the case in most 
of the literature on CMBR polarization. In our case, 
however, the presence of parity-violating sources allows 
nonzero values of the parity-odd correlators and 

and hence these correlators are of special interest 

to us. 

The details of the calculation of and are 

given in Appendix where we find some limiting forms 
of these functions. We have performed numerical eval¬ 
uations for several different values of n and /c* [defined 
in Eq. (^] and assuming g{k) = 1 [defined in Eq. (^]. 
The results for n = 2, n = —2 and n = — 3 are shown 
in Figs. 0-1 The cosmological parameters were set to 
Ha = 0.7, UcDM = 0.25, UBaryons = 0.05 and h = 65 km 
sec“^ Mpc“^. With this choice of parameters, the con¬ 
formal time today is rjQ « 1.5 x 10^ Mpc (we use units in 
which the speed of light and the scale factor at present 
time are set to 1). 

For n > — 1, the angular power spectra of TB and EB 
correlations are dominated by the cutoff. In Fig. |, for 
which we have taken n = 2, the position of the main peak 
in the plots of (l/27r)Z(Z -|- and (l/27r)Z(Z -|- l)Cf^ 

vs I corresponds to the cutoff value of 27r/fc* = 288 Mpc 
or A:*77 o ~ 330. The dotted lines illustrate that spectra 
diverge at large I in the absence of the cutoff. 

As the value of n is decreased, the correlations grad¬ 
ually become cutoff-independent. In Fig. || we plot the 
spectra for n = — 2 with a cutoff at 27r/fc» = 288 Mpc and 
without a cutoff. As can be seen from the plot, above a 
certain scale (Z < 300) the spectra do not depend on the 
cutoff and exhibit a potentially interesting peak struc¬ 
ture. However, spectra still diverge on smaller scales in 
the absence of a cutoff. 

Clearly, any constraint imposed by the future data on 
vq and fc* will be n-dependent. In particular, with the 
chosen form of the initial power spectrum [Eq. ( 0 )] , for 
n > —2 a measurement of or would not con¬ 
strain the helical flow at all. For example, if we optimisti¬ 
cally assume that (l/27r)Z(Z -|- l)Cf^ ~ 10“^°, then for 
n = 2 and 27r/fc* = 288 Mpc we obtain that vq ^ 10^^, 
which is not a very useful bound. For n = — 3 and 
n = —4, with the same value of 27r/A:*, we find vq Y 10“® 
and Vq Y respectively. 

Causality will, in general, constrain the value of the 
spectral index n from below. The bound is obtained by 
setting the real space velocity correlation function to zero 
at causally prohibited separations. This constrains the 
Fourier transform of the correlator to be an analytical 
function of fc, which, in turn, implies that n > 2 for 
small k. Thus, our analysis suggests that the CMBR will 
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FIG. 2: Same spectra as in Fig. for n = —2 with (solid and 
dashed lines) and without (dotted lines) presence of a cutoff 
at 2'Kjkt: = 288 Mpc. 

not be able to constrain primordial kinetic helicity unless 
acausal physics is responsible for producing them. Since 
models that best fit the CMBR temperature power spec¬ 
trum do rely on causality being violated in the past, by 
e.g. inflation or a larger speed of light, it is not incon¬ 
ceivable that the velocity correlations would be acausal 
as well. 


III. KINETIC HELICITY FROM MAGNETIC 
HELICITY? 

In this section we explore the possibility that helical 
magnetic fields may induce kinetic helicity u; at last scat¬ 
tering. As we shall see, the induced velocities are insignif¬ 
icant and can be ignored. 

We shall be interested in the effects of a statistically 
homogeneous and isotropic magnetic field, with possibly 
non-vanishing helicity. If we denote the Fourier ampli¬ 
tudes of the magnetic field by b(k), then, as in Sec. I, 

(6,(k)6,(k')) = (27r)35(3)(k + k') 

X [(5y - kikj)S{k) + ieijikiA{k)\. (17) 

Here S{k) denotes the symmetric part and A{k) the an¬ 
tisymmetric part of the correlator. These functions are 
constrained by [|^ 

5(k) > |A(k)| (18) 

and <^k) > 0, exactly as in the case of kinetic helicity 
[Eq. (^]. Furthermore, 

(B(x) . [V X B(x)]) 


FIG. 3: Same spectra as in Fig. ^ for n = —3. Spectra no 
longer depend on the cutoff scale. 


only depends on A{k) and not on S{k). Therefore A{k) 
represents the helical component of the magnetic field 
and S{k) the non-helical component. 

The Lorentz force F^, due to the magnetic field on the 
(electrically neutral) cosmic plasma is 

FL=jxB, (19) 

where j is the electric current that satisfies Maxwell’s 
equation in the MHD approximation: 

j = V X B . (20) 

Therefore 

(Fi) = (B(x) X [V X B(x)]). (21) 

In the tight-coupling approximation, in which it is as¬ 
sumed that velocities of photons, electrons and protons 
are approximately the same, the Lorentz force induces 
flows of neutral plasma.^ An evaluation shows that the 
quantity (F^) depends on S{k) but has no dependence 
on A{k). In other words, the term with A{k) denotes the 
“force-free” component of the magnetic field. Therefore 
the velocity flow at last scattering is unaffected by the 


^ In the tight-coupling approximation the photons, electrons and 
protons are treated as a single fluid. However, the current j is 
proportional to the difference in velocities of electron and proton 
fluids. It is the Lorentz force due to this slight difference in 
velocities that drives flows in neutral plasma. 
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helical component and the corresponding Doppler signa¬ 
ture on the CMBR can only carry information about S{k) 
and no information about Aik). 

In the preceding discussion, based on the tight- 
coupling approximation, we assumed that the Lorentz 
force acts on an element of the neutral plasma and 
changes its velocity. In reality the coupling of photons 
to electrons is much stronger than that to protons and 
so the plasma at recombination is better treated as com¬ 
posed of two fluids, namely the electron-photon fluid and 
the proton fluid. It is precisely in this approximation that 
the generation of magnetic fields due to cosmic vorticity 
was analyzed by Harrison . 

Our situation is similar to Harrison’s, except for initial 
conditions: we have an initial (force-free) helical mag¬ 
netic field and we need to find the velocity induced by it. 
The analysis, described in Appendix closely follows 
that of Harrison. The result is that the electron-photon 
fluid will gain a vorticity uJe given by 

We = —V^B, (22) 

eue 

where e is the electron charge and Ue is the electron num¬ 
ber density. If we estimate |V^B| ~ B/L?, where L is 
the coherence scale of the field, we find 

|v|-|r-.l-io-'(^)(igi). (23) 

where Bq and Lq denote the magnetic field strength and 
coherence scale at the present epoch and the cosmic elec¬ 
tron number density is ~ 10“®/cm^. (Current bounds on 
cosmic magnetic fields constrain the field strength to be 
less than ^ 10“® G.) Compared to the velocities induced 
by gravitational perturbations (~ 10“^) the velocities in¬ 
duced by helical fields are insignificant. 


IV. A STRATEGY TO DETECT MAGNETIC 
HELICITY 


In the previous section we have shown that only the 
non-helical component of the magnetic field can have a 
signature in the Doppler contribution to the CMBR. If we 
could find another observable that is sensitive to both the 
non-helical and the helical components, we could com¬ 
bine observations and extract the helical component of 
the magnetic field. An observable that does depend on 
both helical and non-helical components is the Faraday 
rotation of linearly polarized sources due to light propa¬ 
gation through a magnetized plasma. 

The CMBR is expected to be linearly polarized and so 
any intervening magnetic fields will rotate the polariza¬ 
tion vector at a rate given by: 

p3 

d9 = X^——aneB-dl , (24) 

STT^mi 


where A is the wavelength of light, a is the scale factor 
normalized so that Otoday = 1, rig is the number den¬ 
sity of free electrons, dl is the comoving length element 
along the photon trajectory from the source to the ob¬ 
server and we are using natural units with h = c = 1 and 
a = e^j(47r) « I/I37. Using the known expression for 
Thomson scattering cross-section. 


fJj' 


Zml 


(25) 


and integrating along the line of sight, we obtain from 
Eq. (H): 


e = 


2'Ke 


f(x) B • dl 


(26) 


where t ( x ) = UeaTCi is the differential optical depth 
along the line of sight, Aq is the observed wavelength of 
the radiation and B = Ba^ is the “comoving” magnetic 
field. 

Faraday rotation depends on the free electron density, 
which becomes negligible towards the end of recombi¬ 
nation. Therefore, the bulk of the rotation is produced 
during a relatively brief period of time when the electron 
density is sufficiently low for polarization to be produced 
and yet sufficiently high for the Faraday rotation to oc¬ 
cur. The average Faraday rotation (in radians) between 
Thomson scatterings due to a tangled magnetic field was 
calculated in Ref. |l^ and is given by 


27re i^Q 


0.08 


Bn 


10-9G 


/ 30GHz 

V 


2 


(27) 


where Bq is the current amplitude of the field and r'o is 
the radiation frequency observed today. 

The amplitude of the CMB polarization fluctuations is 
expected to be of order 10“®, an order of magnitude lower 
than that of the temperature fluctuations. As discussed 
in Ref. Q, detecting a Faraday rotation of order 1° will 
require a measurement which is superior in sensitivity by 
another factor of 10^. Such accuracy is at the limit of 
current experimental proposals but there is a hope that 
it will eventually be accomplished. 

It is usual to define the wavelength independent “ro¬ 
tation measure” (RM) as: 


RM = - / t ( x ) B • dl 

27re ./ 


(28) 


A polarization map of the CMBR at several wavelengths 
will (in principle) give d as a function of A along different 
directions in the sky. From this information the rotation 
measure in any given direction in the sky can be deter¬ 
mined. Hence the polarization map of the CMBR will 
also lead to a “rotation measure map”. Then we can find 
correlations of the RM: 


RR' = (RM(h)RM(n')) , (29) 
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where RM(h) is the rotation measure of the CMBR, ob¬ 
served along the direction ri. Using Eq. dH), this yields 

( 2 ^) 7 ™ 

where 

a = Ki.K;-(k.Ki)(k.K;), (31) 

/3 = ik- (Ki X K;), (32) 

Ki = nJ drj f(7?)e-**^"’', (33) 

K*2 = h' J dr) T{r])e+^^-^'^. (34) 

Eqs. ( ^ ) and (Q) are obtained under the assumption 
that effects of inhomogeneities in free electron density 
along different directions on the sky are of the next order 
in perturbation theory. That allows us to write t ( x ) = 
t(? 7 ) and dl = ndr). 

A crucial feature of RR' is that it depends on both the 
helical and non-helical spectral functions S{k) and A{k). 

The CMBR polarization spectra induced by tangled 
magnetic fields have already been calculated by Seshadri 
and Subramanian |]l^ . They computed the correlator 
between the B-type polarization of the CMBR photons 
commonly denoted by . Defining 

AT^^{l)=Tosl^^^^Cp^, 

where Tq is the CMBR temperature, they found (Eq. (8) 
in Ref. Q) that for I > 1500, 


~ (0.93^K)/| 


(') 

-^rms 

2 

i 

\rJ 

3.10-9G 


1500 


-1-1/2 


(35) 

Here Hrms is the rms magnetic field strength at the 
present epoch and i?* = 770 — is the conformal 
time interval from last scattering to the present epoch. 
The spectral functions S{k) are contained in the mode¬ 
coupling integral I{l/RP as follows: 


X (1-/7^) 


hiq)hi\k + q\)k^ 
m , ^1.1 - 


1 


(k + 2q)i){k + qfj.) 


|k -f- q|‘ 

where |k -|- q| = (fc^ -I- -I- 2 fcg/ 7 )^/^ and 


(36) 


h{k) 


k^Sjk) 


(37) 


The exact form of the Seshadri and Subramanian’s re¬ 
sult is not important for describing our strategy to iso¬ 
late the helical component of the magnetic field. We 
need only observe the crucial feature that AT^^{1), and 
hence CP^, depends on the non-helical spectral function 
S{k) and is independent of the helical spectral function 
A{k). This is because, as discussed in the previous sec¬ 
tion, A{k) is the force-free component of the magnetic 
field and does not induce any velocity in the last scatter¬ 
ing surface. However, it is worth noting a few assump¬ 
tions that enter the analysis in Ref. |^. The first is an 
assumption of Gaussianity by which 4-point functions of 
the magnetic field can be factored into a product of 2 
point functions. The second is the assumption that the 
effects of Earaday rotation at last scattering are not im¬ 
portant. The latter assumption is justified at the high 
frequencies of observation planned for the Planck satel¬ 
lite |l|. 

Hence, if we could use CP^ to obtain S{k) - which 
would only be possible assuming some functional form 
(such as a power law) for S{k) since S{k) occurs within 
some integrals in Eq. ( |^ ) - we could insert the result in 
the expression for RR' given in Eq. (|3^). This will isolate 
A{k) in Eq. (^) and, with some assumptions about the 
functional form of A{k), the cosmic magnetic helicity can 
be evaluated. Even if the assumptions in the derivation 
of Eq. p3) are not completely valid, the RM correlator 
in Eq. ( |30| ) and CP^ have different dependencies on S{k) 
and A{k), and hence the two observations can be used to 
disentangle these two spectral functions. 


V. CONCLUSIONS 


We have analyzed certain P- and CP-violating signa¬ 
tures in the CMBR. If there is kinetic helicity at last 
scattering, it would imprint a signature in the cross¬ 
correlators Cp^ and Cp^. Kinetic helicity can be in¬ 
duced by helical magnetic fields but the effect is too small 
to be significant since the helical component of magnetic 
fields is force-free. Instead we have proposed another 
strategy for detecting the helicity of primordial magnetic 
fields using polarization and rotation measure maps of 
the CMBR. 

The helical magnetic fields produced during elec- 
troweak baryogenesis (~ 10“^^ G at last scattering) are 
several orders weaker than current upper bounds on the 
magnetic field strength (^ I0“® G) at last scattering 

. Therefore the detection of electroweak fields does not 
seem feasible with forthcoming experiments. However, it 
is conceivable that stronger helical fields were produced 
due to some other mechanism and so it is still important 
to think of strategies for detecting primordial magnetic 
fields and helicity. 
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APPENDIX A 

Here we outline our calculation of the signature of ki¬ 
netic helicity in the CMBR. 

In the “total angular momentum” formalism of Hu and 
White 1^, the observable CMB anisotropies T, Q, and U 
are expanded into a set of components 
which are functions of k and rj; here k is the wavenum¬ 
ber in Fourier space and 77 is the conformal time. These 
quantities are such that the power spectra are ex¬ 

pressed through the Z-th components only (see Eq. (56) 
of H), e.g. 


+ (Hvo - v)) 


(A3) 


2Z + 1 = """ 

X (fc,?7)e[=^^^ (fc(?7o - 77)), (A4) 


B 


(± 1 ) 


{Vo,k) 


21 + 1 


I wu 

= —/ dr]fe~'^ 

Jo 

X {k,ri)pl'^^'’{k{r]o-r])),{A5) 


where V is the vector component of the metric per¬ 
turbations, f is the differential optical depth, T{r]) = 
X?" r[i)dr]', and 




(A 6 ) 


The functions satisfy = 


■i.-i -21 


p.-i ,(i) 


Ji' GJP=Jr G and/?)"'=-/3) where 

_ 11(1 + 1) ji (x) 


?(i) 


,(-i) 


3i (x) ^ 


(A7) 


cr = 


dkk 


2 

E 

m— — 2 


(fc,77o) 


(fc,77o) 


21-^1 


21 + 1 


(Al) 


2 


{k,rjo) 


Bt'> (fc,77o) 


= - / dkk^ . , o, , > 

‘ TT J ^ 21 + 1 21 + 1 

m— — 2 

(A2) 

where to = 0 , m = ±1 and to = ±2 denote scalar, 
vector and tensor contributions respectively, and rjo is 
the conformal time today. The quantities Eji^\ 

Pj , in turn, are found from the linearized Einstein and 
Boltzmann equations and are expressed as integrals along 
the line-of-sight of primordial perturbation sources (ini¬ 
tial density and velocity perturbations). Assuming that 
the perturbation sources are random fields with known 
correlations in Fourier space, one can find the expecta¬ 
tion value of any quadratic combination of T, Q and 
U or, equivalently, power spectra (jJ'X',ee,bb cross- 
correlators . 

Scalar perturbations do not generate P-type polariza¬ 
tion. We will also assume that there are no parity vio¬ 
lating tensor sources. The integral solution of the Boltz¬ 
mann equations for the components 0 ^ eI^^\ 

is given by Eqs. (79) and (77) of Ref. 




2Z -h 1 

rvo 

= drj e~'‘ 
Jo 


(^,(± 1 ) {k{r]o-r])) 


jf ix)^ 


+ 1) d f ji (x) 


dx \ x 


_(i) _ \/(Z - 1) (Z + 2) f j[ (a;) ^ ji (x) 


s] (x) = 


+ 


(A 8 ) 

(A9) 


(x) 


V^(Z-1)(Z + 2)j7 ( x ) 

2 X 


(AlO) 


and ji (x) is the spherical Bessel function. It is easy to 
see that, for instance, the correlator depends on the 
parity-odd combination pG)pG) — p(-i)p(-i). 

The GMBR anisotropies sourced by velocity flows are 
predominantly due to the Doppler effect. Since the net 
effect is expected to be small, it is a good approximation 
to keep only terms that are of lowest non-trivial order 
in and in the tight coupling approximation. Then 

Eq. ([A^) becomes 


iVo,k) ^ 
21 + 1 

X 



v^^^Hk,v)jl’^\k{r]o - 77)). (All) 


One can express functions in terms of using 

the linearized Boltzmann equations for temperature and 
polarization anisotropies written in form of infinite recur¬ 
sive series in multipole index Z (see equations (60), (63) 
and (64) of [||). The relevant equations are: 

0p^^ = -Zc^ 0 ^=^^^ - 7-(0p^^ - - V ), (A12) 


























0 


= k 


1 g(±i) _ « 

Vs 7 ^ 


e: 


(± 1 ) 


= -k 


V40 , 


_ TTI' 

7V9 ^ 




^( 0 (±i)p(±i)), 

(A13) 

+ V6P^^^^). (A14) 


If vector metric perturbations are exclusively due to ve¬ 
locity flows, then V is of second order in and its 
effect on CMBR photons is small compared to the in¬ 
duced Doppler shifts. Tight coupling implies that the 
mean free path of photons is negligible compared to the 
scales under consideration: r“^ —> 0 and fc/r <C 1. From 


Eq. (A12), to 0-th order in the tight coupling approxima¬ 
tion, we o btain 0^^^ = To the same order, from 

Eqs. (Add) and ([A^), we find that 


E. 


:±i) 


4 


(A15) 


The quadrupole moments 02 ^^^ and vanish when 

T~^ 0. However, a non-zero dipole term generates 

the quadrup ole mome nts in the n ext lo west order in k/r. 
Using Eqs. ( [A^ ) and ( A15 ), Eq. (Add) can be re-written 
as: 


0 


(± 1 ) 


lv.ey> 


4fc 2 V 2 (± 1 ) _ 

- t/Q - 




3r 7 sVSr 

(A16) 

In the above, the RHS is the source, namely, each of the 
terms on the LHS would be zero in the absence of 0^^\ 
On the LHS, the leading order (in tight coupling) term 


is 0 


(± 1 ) 


Therefore, 0 




p(±i) 


s 4A;r>*^=*=^V(3'\/3'r) and 
3V3r ’ 


(A17) 


From Eqs. (M), (A5), S) and (|A17|) we obtain 
the following expressions for the spectral correlators 


and 

= - ^ /o°° k'^dk k dr]i (AI 8 ) 

_ rii))jl^ {k (ryo - m)) 

X (fc, ryi) (fc, 772 ) — {k, ryi) {k, 772 )) . 


cr = ^ Jo k^dkk^ r dvi r ^^2 

(^k{r]o - rji)) {k{r]Q - 772 )) 

X {k, 771 ) {k, 772 ) - 77^“^) (fc, 771 ) {k, 772 )). 

If we assume that the time-evolution of each Fourier 
mode of the velocity field is independent of k, i. e. the evo¬ 
lution equations contain only k = |k|, then we can write 


the unequal time correlator in Eqs. (AI 8 ) and (A19) as a 
product of the initial power spectrum and the evolution 
functions T{k,ri): 

(fc, 771 ) 7 ;^^^ (fc, 772 ) - f (fc, 771 ) 7 ;^“^^ (fc, 772 )) 

= 7r{k,7ji)T{k,r]2) 

X (k) — {k)v^~^^k)), (A20) 


where (k) — {k)v^~^\k)) is the spectrum 

evaluated at some initial time which we choose to be the 
time of recombination. Using Eq. for the velocity 
power spectrum together with Eqs. ( AI 8 - A20 ), we can 
write Cf ^ and as 

9 /•oo in 

= k^dkvl^g{k)Mk)Mk), (A21) 


9 /•oo in 

C^B ^ _ I k^dkvl-^g{k)Ei{k)Bi{k), (A22) 


TT 


0 


where 


Ai(fc) = / dge '^f jl^{k{go-v))T{k,g) 

JO 


li} + 1 ) r° , jiikjgo-r])) 

2 do ^ k{r]o - v) 

xT(k,ri)fe~^ , (A23) 


/O /*H 0 

Ei{k) = —— dge-^ e^i^\k{go-v))T{k,g) 

yjl- 1)(^ + 1) r° , jijkim-v)) 

3V2 Jo ^ k{go - v) 

x^(e-"r(fc,77)) , (A24) 


fo 

Bi{k) = —— / dge~'^ / 3 j\k{go-g))T{k,g) 


^ j7(fc(77o-77)) 

3V2 Jo ^ Vo-V 

xT{k,g)e-^. (A25) 

As a first step, let us assume that the velocity field is 
simply being red-shifted by the expansion of the universe. 
Then 




(A26) 


For small I we can evaluate and CJ^b approxi¬ 

mately, assuming g{k) = 1 : 


A,(7-) » (A27) 


Hk) 


Vi^- 0(1 + 1 ) jlikgo) 


3 V 2 


kgo 


1 + 2 ( — 
ar 


Bfi\^ V(^- l)(^ +1) ji(fc??o)^-i 

B>i\K) ~ ^—7=- 


3 V 2 




(A28) 

(A29) 
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cT^ = 


(Z + l)v/?(I^ vl 

StT '^recVO 

nOO 

/ dxx'^+'^jfix) , 


{k*rio) 


-(ra+3) 


(A30) 




(Z-1)(Z + 1) 1^0 ^ ^_(„+3) 


9tt 


1 + 2 1 — 

ar 


'^rec'^0 


-{KmY 


dxx'^^ jf(a;).(A31) 


The last integral can be evaluated: 

r, »+.- 2 . > 2"-‘ir(-n)r(i + i + f) 

I -i” 2.W= + 

We find for n < 0 that 


. (A32) 


CTB,EB ^ p+2^ 


(A33) 


For n > 0 the integral diverges and is dominated by the 
range near the cutoff: 



dxx^+^jf{x) 


n+2 


n + 1 


(A34) 


TB.EB 


and C'l 

The expressions in Eqs. ( A21 ) and ( |A22| ) have been 
evaluated numerically and the results are described in 
Sec. g. 


APPENDIX B 


Here we find the vorticity induced by magnetic fields at 
last scattering in the two-fluid approximation where the 
photon and electron velocities are equal, but the proton 
velocities can be different. This is closely analogous to 
Harrison’s calculation of the induced magnetic field due 
to cosmic vorticity [|ll|| . 

The starting point is the Euler flow equations for the 
velocities v;,, Vg, v.y of the baryon {i.e. proton), electron 
and photon fluids. These equations are all of the form 



EPb 


where 


dv 

dt 


dv 

'm 


+ (v • V)v 


is the total derivative of the velocity. 


Fi = n^Zie(E + x B) 


(HI) 


(B2) 


(B3) 


is the force from the electromagnetic field per unit vol¬ 
ume, ZiC is the charge of particles constituting the fluid 
i, and = —Pji is the rate of momentum transfer to 
the fluid i due to collisions with the fluid j [hence the 


negative sign in Eq. (Bl)] per unit volume. For baryons 
we take Zf, = +1. We can also add a gradient of the 
gravitational potential to the equation for baryons. 
The density of particles is pi = riimi. 

Collisions between photons and baryons are disre¬ 
garded because they transfer much less momentum than 
collisions between other fluids. Collisions of electrons 
with baryons determine the conductivity a of the plasma: 


erie. 

Peb = - J- 

a 


(B4) 


The velocities of the fluids are almost equal to each 
other due to tight coupling, however the velocities of elec¬ 
trons and photons are closer to each other than the ve¬ 
locities of electrons and baryons: 


|v..y - Vej < jVe - Vbj < Vfe. (B5) 


We take the velocities of the electrons and photons to be 
the same, except when computing Pe-y- 

We now rewrite the Euler flow equations for baryons 
electrons, expressing the combination E + v^ x B through 
other quantities: 

E + VfcxB = ^(^+V(^) +ij, (B6) 
e \ dt J ^ 


„ We dvp 1. 

E + Ve X B = - — H— J 

e dt a 

The flow equation for photons is 


67 • 


(B7) 


(Pt +P 7 ) ^ “ (P7 +P7)^</>- Pe7: (B 8 ) 

and we use the equation of state = p^/d. 

Maxwell’s equations in the MHD approximation are 


/9R 

VxE = - —, Vxj = -V2B. 
at 


(B9) 


Taking curl of Eq. (^), using V • B = 0 and Eq. (|B2| ) 
for the total derivative d/dt, we obtain 


e e 

-V^B = — 

mb 


mbcr 


dB 

dt 


■ B (V • Vb) - (B • V) Vb 
duJb 


+V X [(vb • V)vb] - (vb • V) u)b + —. (BIO) 


Next we write 


Vb = u + w. 


(Bll) 


where u = Rr/R is the velocity due to Hubble expansion, 
R{t) is the universal scale factor, and w is the vortical 
velocity in the angular direction around the line of vor¬ 
ticity which is taken to coincide with the direction of the 
magnetic field lines (assumed to be parallel). Then, after 
some manipulations, we obtain 



2B- = -V^B. (B12) 

R a 
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We can rewrite this as 


^ + brA = —v^b. 

at \ e J a 


(B13) 


In Eq. (B7), Harrison neglects the dve/dt terms and 
also takes Ve = in the Ve x B term. Taking curl of 
that equation then gives 

dB R I 1 , , , 

—+ 2B- =- VxPe^ + -V2B. B14 

at R eUe cr 

(Harrison also neglects the V^B term which we ke ep 
here.) The unknown quantity V x Pg-y in Eq. (B14) is 
found by taking the curl of Eq. (B8), 






R 




= —3V X Pe 


(B15) 


Combining Eqs. (B14) and ( B15| ), we get 

dB ^R 4 1 f du}^ r\ 1 
dt i? 3 erip ^ \ dt ^ R I o 


V^B. (B16) 


Note that p ^t) = p^^R~'^ and Pe (t) = pi^^R~^. Multi¬ 
plying Eq. (B16) by R^ and rewriting it as a total time 
derivative, 


d 

dt 


BR^ - 


,2 


?>ep, 


(0) 


i?2 

= —V^B. 

cr 


(B17) 


Subtracting Eq. ( [BITD from Eq. ( |B13D and integrating 
with the initial condition of zero vorticity finally gives: 


uibR^ Ap^^meR 

- UJl) — - 777 ^- it). 


( 0 ) 


Sep, 


Additionally, we take the curl of 


Vfe - Ve = -j 

erie 


equate u:-y and We, and obtain 


LiJb — eiJ-y = -V^B. 

eUe 


7 • 


(B18) 


(B19) 


(B20) 


No w we have two simultaneous equations, Eq. (BIS) 
and ( B20 ) for the two vorticities Wg = u)-y and Wf,. Rec¬ 
ognizing that Pj <C Pb at last scattering, we then solve 
the simultaneous equations to get 


U3e = -V"^B . 

erie 


(B21) 
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